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With increasing electron-phonon coupling as described within the r-J-Holstein model, bipolaron 
kinetic energy is lowered in comparison with that of the polaron. This effect is accompanied with 
"undressing" of bipolaron from lattice degrees of freedom. Consequently, the effective bipolaron mass 
becomes smaller than the polaron mass. Magnetic as well as lattice degrees of freedom cooperatively 
contribute to formation of spin-lattice bipolarons. 
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I. INTRODUCTION 

In strongly correlated systems superconductivity (SC) 
may occur as a consequence of the kinetic energy lower- 
ing as opposed to a standard BCS-type superconductor 
where there is a slight kinetic energy raise accompanied 
by the lowering of the potential energy that stabilizes 
the SC state^™. Experimentally, kinetic energy lowering 
would reflect in the violation of the low-frequency optical 
sum rule and lead to a change of high-frequency optical 
absorption upon entering SC stated 

The idea of in-plain kinetic energy driven pairing, orig- 
inally proposed by J.E. Hirsch, is based on the Hubbard- 
type model where hopping depends upon the occupation 
numbe r 2 . The possibility of kinetic energy gain has been 
investigated as well in more general models with corre- 
lated electrons, such as the standard Hubbard model. 
Using variational Monte Carlo method authors of Ref. 6 
find kinetic energy gain in SC state with a d-wave symme- 
try above a critical value of U c . More recent calculations 
based on the dynamical cluster approximation show that 
pairing is driven by the kinetic energy gairP. 

The main motive for the kinetic energy driven pairing 
in models possessing at least short range antiferromag- 
netic correlations relics on the argument that the mo- 
tion of the single hole is obstructed due to formation of 
strings of misaligned spins left in the wake of the prop- 
agating hole. A pair of holes that propagates coherently 
should lower its kinetic energy as one hole moves in the 
wake created by the other hole. This naive argument 
was challenged by S.A. Trugman® who suggested, that 
a pair of holes is less mobile than originally anticipated 
when assuming the simple string argument. Lower pair 
mobility occurs due to a frustration effect, which arises 
from the fermion exchange processes. Cluster dynamical 
mean field studies^- on the t-J model nevertheless show 
a small kinetic energy gain in the underdoped regime as 
the system enters SC state, while there is a slight kinetic 
energy raise in the overdoped regime. Authors of RefP^ 
have demonstrated the existence of the kinetic energy 
driven superconductivity in the t-J model using the spin 
polaron technique. 



A scenario of the kinetic energy gain upon pair for- 
mation is inherently connected with increased pair mo- 
bility and consequently with lowering of the effective 
mass. In Holstein-type models, however, the effec- 
tive bipolaron mass is in the strong electron-phonon 
(EP) coupling regime typically much larger than the 
polaron one. A heavy bipolaron effective mass repre- 
sents one of the main obstacles for bipolaronic theory of 
superconductivity^] Addition of Coulomb interaction^ 
and generalization to physically more relevant Frohlich 
type EP interactional!^ contribute to a substantial de- 
crease of the bipolaron mass. Nevertheless, the bipolaron 
remains heavier in comparison with the polaron. 

While it is widely accepted that strong correlations 
govern the physics of high-T c superconductors^, the no- 
tion of the importance of lattice effects with the empha- 
sis on their role in formation of the SC state is as well 
gaining momentum^. I n this letter we compare physical 
properties of systems with one and two holes coupled to 
quantum phonons doped in the Heisenberg antifferomag- 
net. 

We first show that there is no kinetic energy gain upon 
bipolaron formation in the pure t-J model. Throughout 
this paper we use the term kinetic energy for the expec- 
tation value of the hopping term in the t-J model. We 
should point out that this term represents only the ki- 
netic energy of the lower Hubbard band, as already noted 
by Wrobel et al. in Ref. 10 . Switching on EP coupling we 
discover that the interplay between the kinetic, magnetic, 
and elastic energy leads to a formation of a spin-lattice 
bipolaron that in the crossover regime between the weak 
and strong EP coupling gains the kinetic energy relative 
to its polaron constituents. In this regime the effective 
mass of the spin-lattice bipolaron is decreased relative to 
the effective mass of the polaron. 



II. MODEL AND METHOD 

We employ a recently developed method to solve sys- 
tem of one and two holes in the t-J model defined on 
an infinite two-dimensional lattice, coupled to lattice de- 



2 



grees of freedorrflSHSI] We investigate the influence of 
EP coupling via the simplest extension of the t-J model 
where holes couple to dispersioneless phonons. The 
model describes the influence of apex oxygen vibration 
on doped holes, propagating in lightly doped CuO planes 
in cuprates: 



H = -* E (4a*+ h - c -) + J E( s i s J - \ n ^) 
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where Ci jS = Ci. s (l — ni._ s ) is a projected fermion op- 
erator, t represents nearest neighbor overlap integral, 
the sum runs over pairs of nearest neighbors, ai 

are phonon annihilation operators and rii — 2W n »i«' g 
and ujq represent EP coupling constant and the Einstein 
phonon frequency, respectively. 

While the numerical method has been in detail de- 
scribed in previous workfl^ID we om y briefly highlight 
a few most relevant elements of the method. The con- 
struction of the functional space for one and two holes 
starts from a Neel state with one or two holes located 
on neighboring Cu sites and with zero phonon quanta. 
In the case of a high symmetry point at k = (0,0), the 
parent state of two holes can be chosen to exhibit a point- 
symmetry, belonging to a particular irreducible represen- 
tation of the point group C^ v . Such a state is expressed 
as |0(°.°)) a = ^ 7 (-l) M »(T)c c T |Neel;0), where sum runs 
over four nearest neighbors in the case of d— and s— wave 
symmetry and over two in the case of p^j^-wave while 
M a (7), a 6 {d, s,p} sets the appropriate sign. 

We generate new parent states by applying the gener- 

ator of states {\^ NhM) ) a } = (H kin + H?) Nh \^) a 

where -Hkin represent the first term in Eq. [I] H g rep- 
resents third term in Eq. [I] When parameter M > I 
is chosen, the functional generator creates states with 
additional phonon quanta. This approach ensures good 
convergence in the strong EP coupling regime where the 
ground state contains multiple phonon excitations^. In 
most cases we have used Nh = 6 and M = 8 that lead 
to N st = 24 x 10 6 states. Full Hamiltonian in Eq. [l] 
is diagonalized within this limited functional space tak- 
ing explicitly into account translational symmetry. As of 
now we refer to one and two hole states as polaron and 
bipolaron, where polaron (bipolaron) signifies a hole (two 
holes), dressed with spin as well as lattice excitations. 



III. RESULTS 
A. t-J model 

We first investigate the possibility whether within the 
framework of the pure t-J model a bipolaron gains the ki- 
netic energy in comparison with a polaron. To this effect 
we show in Fig. [lja) the expectation value of the kinetic 
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Figure 1: (Color online) (a) (i/kin) / Nhoi for the AW = 1 (po- 
laron) and 2 hole (bipolaron) in the t-J model with A = 0. 
There are three nearly overlapping curves for each system ob- 
tained using the following parameters Nh — 10, 12, and 14 
with M = in all cases that lead to the following numbers of 
basis states: N st = 3.1 x 10 5 ,2.6 x 10 6 , and 21 x 10 6 for the 
TVhoi = 2 case. Numbers of polaron states (AW = 1) are typ- 
ically a factor of 4 smaller; (b) Akin vs. J/t for the t-J model 
(full curve) and for the anisotropic t-J z model (dashed line). 
A discontinuity in the latter case marks the crossover from 
the p— wave bipolaron ground state at small J z /t < 1.2 to 
d— wave at large J z /t; (c) comparison of (ii/kin}/AW between 
AW = 1 and 2 systems vs. A. In this and all subsequent 
plots we used generator of states with Nh = 6 and M — 8 
that led to a Hilbert space with AW = 24 x 10 6 states and a 
a maximal number of phonons AW = 48, (d) Aki n represent- 
ing the difference of kinetic energies between bipolaron and 
polaron as defined in the text vs. A for different values of the 
exchange interaction J/t. The thicker curve is for J/t = 0.4. 



energy per hole, i.e. (Hun)/N\ i oh vs - J/t °f the one and 
two-hole Hilbert space. Except in the unphysically small 
J/t < 0.15 we find no kinetic energy gain of the bipo- 
laron state. This result is in agreement with Trugman's 
suggestion^. The frustration that arises from the fermion 
exchange processes impedes bipolaron motion. There is 
no such effect in the polaron case. 

A higher mobility of the bipolaron should be more pro- 
nounced in the case of the t-J z model due to a lack of 
spin-flip processes that erase pairs of overturned spins 
and contribute to the gain of the kinetic energy of po- 
larons in the isotropic case. For this reason we have 
investigated the difference between kinetic energies per 
hole between bipolaron and polaron Akin = (if kin) ^ /2 — 
(iJkin)^ 1 '' in a wider, even though unphysical range of J/t, 
see Fig. [ijb). We find A kin < for J z /t > 1.5. The dis- 
continuity is a consequence of a crossing between p— wave 
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symmetry of the pair at small J z /t to d— wave symmetry 
for J z /t > 1.2. In contrast, in the isotropic t-J model 
Akin > in the whole expanded J/t regime, presented in 
Fig. Qb). 

On the more technical side we report on a test of 
the convergence of our method. In Fig. [lja) we present 
nearly overlapping curves obtained using three different 
Hilbert spaces, generated by — 10, 12, 14 and M = 0. 
Note, that in our calculation the maximal allowed hole 
distance is: l max — Nh + 1 = 15 in the case of the largest 
Hilbert space with no phonons. This should be compared 
with exact diagonalization calculations on finite square 
lattices with N sites, where l max = \J N/2 = 4 in the 
case of N = 32 sitesPD. 



B. t-J-Holstein model 

We now switch on the EP coupling. In Fig.[TJc) we plot 
the kinetic energy per hole vs. dimensionless EP coupling 
constant A = g 2 /8u>ot. At small A kinetic energies of 
polaron (iVhoio = 1) and bipolaron (iVhoie = 2) increase 
linearly with A. Such behavior is characteristic for the 
weak EP coupling regime. In the regime 0.22 < A < 0.32 
the kinetic energy (per hole) of bipolaron crosses below 
the kinetic energy of polaron. The onset of this regime 
coincides with the crosso ver t o strong coupling regime 
of the spin-lattice polarorP^U. If we roughly define the 
crossover to the strong EP coupling regime as a point of 
the steepest increase of (-Hkin) vs. A, we discover, that the 
polaron state enters strong EP coupling regime at smaller 

A^ ~ 0.22 than the bipolaron state, where ~ 0.28. 
The two inflection points are as well indicated with open 
circles in Fig. [ljc) . 

The dependence of the kinetic energy gain on the mag- 
netic exchange interaction J/t is investigated in Fig.[l|d) 
where we follow Akin vs. A using different values of the 
exchange interaction. As J/t increases the kinetic en- 
ergy gain is reduced. This behavior is in contrast with 
the anisotropic case where at A = 0, Akin becomes nega- 
tive only at large J z /t, see Fig. [TJb). The kinetic energy 
gain, observed in the isotropic case, is driven by the EP 
coupling. It is positioned in the crossover from weak to 
strong EP coupling regime. In terms of the magnetic ex- 
change interaction it is located well within the physically 
relevant regime of J/t € [0.3,0.4]. 

In Fig. [2ja) we present the average hole distance 
(d) = J2 r r -P( r )' OI a bipolaron state, where P(r) repre- 
sents the probability of finding a hole-pair at a distance 
of r: P(r) = (£ <i?y) n*r$5 [|i - j| - r])/<E<i*j> nfnf). 
At J/t — 0.4 and A = two holes form a bound bipo- 
laron with the largest probability at a d istance o f r = v2, 
as consistent with previous calculation d 21 * 23 * 2 -^. With in- 
creasing EP coupling A, (d) experiences the steepest de- 
crease in the middle of the regime where Akin < 0, i.e. 
for 0.22 < A < 0.32. This is reflected in the change 
of the shape of bipolaron where the hole distance with 
the largest probability P(r) crosses over from r = \/2 at 



I ' I 1 I ' I 

| O — O k=0.00 

1X=0.22 - 
>k=0.32 1 




0.1 0.2 0.3 X 0.1 0.2 0.3 I 



Figure 2: (Color online) (a) Average hole distance vs. (d). 
Vertical dashed line indicates the onset of the regime of the 
kinetic energy gain, Akin < as seen in Fig. [jja). (b) Prob- 
ability P(r) of finding holes at a distance of r computed at 
three different values of A = 0, 0.22 and 0.32, the latter two 
coinciding with positions of two arrows in (a), (c) (n p h)/-Nh i 
vs. A of polaron (full line) and bipolaron system (dashed 
line). The dot-dashed line represents the fit of the polaron 
result to a straight line in the regime A > 0.375 that is given 
by y = C + 8.8\t/wo (d) effective masses per hole and 
m' 2 - 1 of polaron (full line) and bipolaron (dashed line) system, 
(e) and (f) R n and _R m ratios between phonon numbers and 
effective masses as presented in (c) and (d) respectively (see 
also definitions in the text). 



A = to r = 1 at A = 0.32, see Fig. |2jb). It is somewhat 
counter-intuitive that such shrinking of the bipolaron size 
simultaneously gives rise to Akin < 0. Close proximity of 
holes namely restricts the range of hopping that conse- 
quently leads to a raise of the kinetic energy, as as well 
seen in Fig. [T|c). 

We gain additional insight into the mechanism lead- 
ing to Aki n < by presenting comparison of the average 
number of phonons per hole between polaron and bipo- 
laron. In Fig. [2tc) we plot (n p h)/Ahoi vs. A. At small 
A the average (n p h) /A no i is for bipolaron slightly larger 
from that of the polaron, while at A ~ 0.22 it crosses be- 
low the polaron result. In the large A limit both expecta- 
tion values should approach the strong coupling regime. 
Indeed, (n p h) /A no i for polaron for A > 0.27 approaches a 
straight line characteristic for the strong coupling result: 
(« P h) = g 2 /^o = 8A</w , see the fit in Fig. [2jc). The 
same is not true for the bipolaron, where (riph)/N\ lo \ is 
below polaron result for A > 0.22 and finally approaches 
the strong coupling result above A > 0.35. We should also 
note that in the Akin < regime, say around A = 0.26, 
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(^ph) /iVhci for the polaron system exceeds result for bipo- 
laron by nearly 50%. Results are consistent with the ob- 
servation that the bipolaron crosses over to the strong 
coupling regime within a wider crossover regime and at 
larger A than the polaron. 

In connection with Fig. [2^c) we note that in the vicin- 
ity of the physically relevant values of A 6 [0.22,0.32] 
the average number of phonons in the bipolaron state is 
( n ph) < 20. This is to be compared with the maximum 
number of phonon quanta contained in the Hilbert space 
N ph = N h * M = 48. We can conclude that Hilbert 
space used in our method contains sufficient amount of 
phonon degrees that empowers our calculation reach- 
ing full convergence. To further investigate the differ- 
ence in the phonon number between bipolaron and po- 
laron state we plot in Fig. [2je) the corresponding ratio 
Rn = (™ph) < - 2 V2(™ph) < ' 1 ' that shows 40% decrease in the 
avergage phonon number per hole in the regime Akin < 0. 
Note that in the strong coupling regime R a —> 1. 

We proceed by presenting comparison of effective po- 
laron and bipolaron masses. The polaron dispersion of 
the t-J model has a minimum at k = (it/2,it/2) and is 
highly anisotropic. We compute the effective mass ten- 
sor in its eigendirections m aa = 2t (9 2 i?(k)/<9k<9k) Qc[ . 
Taking into account anisotropic dispersion we define the 
polaron mass as to/ 1 * 1 — ^/to,||TO_l, were to,|| and m± are 
effective polaron masses along nodal and anti-nodal di- 
rections respectively. In contrast, bipolaron has the en- 
ergy minimum at k — with locally isotropic disper- 
sion. We compute effective bipolaron mass per hole from 
to' 2 * 1 = m xx /2. In Fig. [2^d) we present to/ 1 * 1 and to/ 2 - 1 
vs. A. In the weak coupling regime we find the expected 
result where to' 2 * 1 > to/ 1 '. For A > 0.23 the opposite be- 
comes true. Results in Fig.^d) represent, at least to our 
knowledge, the first example where to/ 2 * 1 < to/ 1 - 1 . The ra- 
tio R m = to' 2 - 1 /to/ 1 - 1 in Fig. j^f) drops down to R m ~ 0.2 
around A ~ 0.3. This result is highly unusual. For com- 
parison we draw attention to a well known result for the 
Holstein model where the effective mass of two particles 
forming a lattice bipolaron singlet, scales in the strong 
coupling regime as to/ 2 ' oc exp[4(g/wo) 2 ] in comparison 
to to/ 1 - 1 oc exp[(g/u;o) 2 ] that leads to to/ 2 * 1 >> to/ 1 - 1 . 
In the Holstein Hubbard model, however, the on-site 
Coulomb interaction gives rise to a formation of an inter- 
site bipolaron with an effective mass that is comparable, 
nevertheless always larger than the polaron mass^l. Re- 
cently Hague et alf^ examined a bipolaron defined on 
a triangular lattice and found unusually small effective 
bipolaron mass due to a crablike motion. Nevertheless, 
bipolaron effective mass remains larger than that of the 
polaron. 

A note of caution: as seen in Fig. [2^d) the absolute 
value of the effective bipolaron mass is rather large all- 
ready at the onset of the Akin < regime, to/ 1 - 1 ~ to/ 2 - 1 ~ 
60 (in units of free electron mass) at A = 0.23. This 
shortcomming can be easily alleviated by taking into 
account physically more relevant longer-range Frohlich 
EP interaction that may reduce to/ 2 -' and m/^up to an 



order of magnitude^^. Experimental results on the 
efective mass seem to be slightly ambiguous. Neverthe- 
less, latest de Haas-van Alphen measurements of the cy- 
clotron effective mass inside of the superconducting dome 
of YBCO show divergence with decreasing doping. They 
have measured effective masses as large as 4.5 free elec- 
tron masses^. 
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Figure 3: (Color online) (a) (//kin) /AW for AW = 1 and 
2 systems vs. A, (b) average hole distance (d) vs. A, (c) 
(nph) /AW vs. A of polaron (full line) and bipolaron system 
(dashed line), (d) effective masses per hole rrv- 1 ' and nv' for 
polaron (full line) and bipolaron (dashed line) system, (e) and 
(f) -Rn and Rm ratios between phonon numbers and effective 
masses as presented in (c) and (d) respectively. Parameters 
of the model in all pictures (a, f) are: J/t — 0.4 and 

LUo/t = 1.0. 

In our search for deeper understanding of the phonon 
driven kinetic energy lowering and bipolaron mass renor- 
malization we have investigated as well the regime with 
larger uj , i.e. uo/t = 1.0. In this case the effect of 
EP coupling on the kinetic energy lowering disappears as 
seen from Fig. |3ja) even though around A ~ 0.4 Akm 
approaches zero, see also the insert of Fig. |3ja). On 
the other hand the effective bipolaron mass still shows 
a slight decrease with respect to the polaron mass in the 
regime 0-35 < A < 0.6, as shown in Figs, -^d) and (e). 
This effect is closely connected to the average number 
of phonons that shows barely detectable decrease in the 
same parameter regime, see Figs. J3jd) and (e). 



C. Attractive t-J-V model 

We have tried to reproduce the kinetic energy lovering 
and the effective bipolaron mass renormalization using 
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Figure 4: (Color online) (a) Akin vs. V/t for the t-J-V 
model, (b) effective masses per hole mS^ and m' 2 - 1 for po- 
laron (dashed line) and bipolaron (full line) system vs. V/t, 
and (c) average hole distance (d) vs. V/t. We have used basis, 
generated by iV n — 12 and M — 0. 

a t-J-V model where the effect of phonons is replaced 
by an effective nearest neighbors attractive interaction 
V. To this effect we have added an attractive term of 
the form Hy = — VJ^/yv nfnj 1 where nj 1 represents the 
hole-density operator, to a standard t-J model. From 
Fig. Qa) we conclude that with increasing V/t Akin re- 
mains positive and monotonically increases. We find no 
kinetic energy gain in this simplified model. Neverthe- 
less, we find a decrease of the effective mass of the bipo- 
laron system as presented in Fig. [4jb) . Effective mass 
ration reaches its minimum value R m ~ 0.75 around 
V/t = 3.0. This should be compared to R m ~ 0.13 
in the case of oj = 0.2 and A = 0.3, Fig. gf). Even 
though the effect of nearest neighbor attraction V on the 
average distance (d) is similar to the effect of increas- 
ing A (compare Figs [2]Ja), [3]^b), and He)), the influence 
of the EP coupling on hole motion in the t-J model can 
not be entirely explained using an effective attraction be- 
tween holes. Using the simplified model we nevertheless 
discover an important mechanism whereby the effective 



bipolaron mass decreases at fixed value of J ft as the aver- 
age distance between holes decreases by invreasing V . In 
contrast to spin lattice bipolaron where its cumulative ef- 
fective mass at u>o = 0.2 and A ~ 0.3 becomes lower than 
that of a polaron, R m ~ 0.13 < 0.5, spin bipolaron with 
attractive hole-hole interaction remains heavier than the 
polaron since R m ~ 0.75 > 0.5. 

IV. DISCUSSION 

The main mechanism behind the gain of the kinetic en- 
ergy as well as lowering of the effective mass of bipolaron 
emerges from the competition between kinetic, magnetic 
and lattice degrees of freedom. As two holes form a bipo- 
laron in the regime where Akin < 0, the system minimizes 
the kinetic energy at the expense of increased elastic en- 
ergy. The gain of the kinetic energy contributes to energy 
splitting between states of differnt point-group symme- 
tries. This is in agreement with recent results of RefPHl 
where it has been shown, that EP coupling to transverse 
modes in-plain oxygen vibrations stabilizes d— wave sym- 
metry. The increase of the elastic energy emerges as un- 
dressing of bipolaron from lattice degrees of freedom as 
seen as the decrease of (n pn )/iVhoi below the respective 
polaron values. The gain in the kinetic energy does not 
represent the "glue" for the formation of the bipolaron 
since for J/t > 0.2 bipolaron is already formed in the pure 
t-J model where no gain in the kinetic energy is found. 
Instead, it emerges as a side product of the competition 
between the kinetic, magnetic and elastic energies. 

Our results lead to a novel paradigm where in a cor- 
related system, coupled to quantum lattice degrees of 
freedom, upon pair formation the bipolaron mobility in- 
creases due to a lower effective mass as well as due to 
a detectable gain in bipolaron kinetic energy. The at- 
tractive potential for binding of bipolaron appears as a 
cooperative interplay between magnetic and lattice de- 
grees of freedom. While the original idea of the kinetic 
energy gain as a mechanism for hole-pairing has been 
proposed for a system of correlated electrons, our finding 
opens the possibility where lowering of the kinetic energy 
in a correlated model is driven by the EP interaction. 
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